Quantum superconducting criticality in graphene and topological insulators 



Bitan Roy 1 , Vladimir Juricic 2 , and Igor F. Hcrbut 31 

11 National High Magnetic Field Laboratory, Florida State University, Tallahassee, Florida 32306, USA 
2 Instituut-Lorentz for Theoretical Physics, Universiteit Leiden, 
P.O. Box 9506, 2300 RA Leiden, The Netherlands 
3 Department of Physics, Simon Fraser University, Burnaby, British Columbia, Canada V5A 1S6 

The field theory of the semimetal-superconductor quantum phase transition for graphene and 
surface states of topological insulators is presented. The Lagrangian possesses the global U(l) sym- 
metry, with the self-interacting complex bosonic order-parameter and the massless Dirac fermions 
coupled through a Yukawa term. The same theory also governs the quantum critical behavior of 
graphene near the transition towards the bond-density-wave (Kekule) insulator. The local U(l) 
gauged version of the theory which describes the quantum semimetal-superconductor transition in 
the ultimate critical regime is also considered. Due to the Yukawa coupling the transitions are found 
to be always continuous, both with and without the fluctuating gauge-field. The critical behavior 
is addressed within the dimensional regularization near four space-time dimensions, and the calcu- 
lation of various universal quantities, including critical exponents and the universal mass-ratio is 
reported. 



Dirac quasiparticlcs represent low-energy excitations 
in various low-dimensional condensed-matter systems, 
such as graphene and topological insulators (TIs). In 
graphene, pseudorelativistic Dirac quasiparticles emerge 
from hopping of the electrons on the underlying hon- 
eycomb lattice? 1 ^ while on a surface of a TI, they re- 
sult from an odd number of band inversions in the bulk 
of the system leading to a topologically non-trivial Z2 
invariant.— In both cases, the linearly dispersing Dirac 
quasiparticles give rise to a semimetallic ground state, 
stable against weak electron - electron interactions. 4 

When the repulsive interactions are sufficiently strong, 
however, a plethora of insulating phases can in principle 
be realized in graphene4~— Furthermore, Dirac fermions 
in graphene can also condense into four different gapped 
superconducting states, if the net interaction acquires an 
attractive component The simplest of them, which 
will be the subject of the present study, is the uni- 
form, spin-singlet s-wave pairing, favored by a sufficiently 
strong on-site attractive interaction.— Two of the re- 
maining pairing gaps are spatially inhomogeneous spin 
triplets, which break the translational symmetry of the 
honeycomb lattice into Kekule patterns^ They are fa- 
vored by the sufficiently strong nearest-neighbor attrac- 
tion. Finally, yet another triplet pairing with an /- 
wave symmetry can be stabilized by a strong second - 
neighbor attraction.— On the other hand, due to their 
reduced number of fermionic components, the massless 
Dirac fermions residing on the surface of TIs can acquire 
a superconducting gap only by pairing into an s-wave su- 
perconducting state. Possible inhomogeneous^ and chi- 
ral d + id superconducting states in doped graphene,- 
as well as proposed realizations of Majorana fermions 
in graphene and TIs,— ~— make the study of supercon- 
ducting instabilities of Dirac fermions in low-dimensional 
condensed-matter systems theoretically and experimen- 
tally interesting and timely! 14 ' 15 

The bosonic order parameters (OPs), characterizing 
both the insulating and the superconducting states, are 



composite objects of Dirac fermions, and may exhibit 
different symmetries. Besides the usual self- interaction, 
the OPs here are also coupled to the massless Dirac 
fermions via Yukawa term^ Previously, we studied the 
Ising and Heiscnbcrg universality classes of the transi- 
tion into charge-density-wave and spin-density-wave in- 
sulators, respectively. The Cooper pairs are of course 
charged, and consequently the corresponding field the- 
ory possesses a global U(l) symmetry. We therefore here 
develop the U{1)- symmetric field-theoretical description 
of the quantum semimetal-superconductor transition in 
graphene and surfaces of TIs. All the coupling constants 
in this effective theory are marginal in d = 4 space- 
time dimensions, which enables us to perform the e— 
expansion with e = 4 — d to address its critical behav- 
ior. We find that the semimetal-superconductor transi- 
tion is always continuous, and compute the critical ex- 
ponents associated with the transition; in particular, the 
correlation-length exponent and the anomalous dimen- 
sions for both the OP and the fermion fields are ob- 
tained. Besides describing the semimetal-superconductor 
quantum phase transition, the U(l) field theory should 
also pertain to the quantum phase transition from the 
semimetallic into an insulating state with the dynami- 
cally generated Kekule mass, which breaks the transla- 
tion symmetry of the latticed Motivated by this physi- 
cal problem, as well as by the theoretical possibility of a 
fully Lorentz invariant semimetal-superconducting tran- 
sition in these Dirac systems, we extend our theory to in- 
clude a fluctuating U(l) gauge field, that would describe 
the coupling of the electromagnetic field to the bosonic 
OP and the Dirac fermions near the critical point The 
Lorentz-symmetric critical point describing the continu- 
ous transition in this theory is charged, and also turns out 
to be stable for any physical number of fermion flavors 
(that is, for N f > 0.142). 

To set up the problem, we first consider the 
pairing of the the gaplcss excitations in graphene, 
around the two incquivalent Dirac points, at 



±K, described by an eight-component Dirac- 
Nambu spinor *t(fc) = [* 1 f (fc),* t _(fc)], where 
tft(fc) = [ui(k),vl(k),au_ a (-k),av_ a (-k)]M Here, 

K = (1, l/v3)(27r/av3), with a being the lattice con- 
stant, k = (w,k) is the three-momentum and k = K + q, 
|q| <C |K|. cr = ± is the spin projection along the z-axis. 
u a and v a are the Grassmanian fields on two sublattices. 
The free Dirac Lagrangian in this representation assumes 
relativistically invariant form L° = i*f?(x)ao <8 7 At (9 At , J , (x), 
where = / d 3 ke ikx f(k), fj, = 0,1,2, x = (r,r) 

with t as the imaginary time and summation over the 
repeated indices assumed. The 7-matrices are defined as 
70 = cr 3 <g> (73, 71 = cr (8 cr 2 , 72 = 00 ® 01, 73 = 01 ® t7 3 , 
and 75 = <T2 <E> CT3, where {00, cr} form the Pauli basis 
for two-dimensional matrices, and we take if? = ^jo, as 
usual. The s-wavc superconducting OP reads 

$(x) = (*t (x)a <8 («7o73 cos(p + i7o7s sin y>)*(ar)), (1) 

with 93 as the phase of the superconducting OP. The OP 
anticommutes with the particle number operator ./V = 
Co ® *7375 and commutes with all the three generators 
of the spin rotations, S = cr (g> 74, and hence represents a 
spin-singlet. Moreover, it is even under the exchange of 
the sublattices, or of the Dirac points. 

In terms of Nambu's (particle-hole doubled) spinor ba- 
sis, vfft = (ct , ct , c r. , — c_r . ), the Lagrangian for the 

gapless surface states of TIs also adopts the relativistic 
form L°j = iif?(x)'^ ll d ll 'f'(x), with the 7-matrices defined 
as 70 = cr 3 <8> 0-3,71 = ® cr 1 ,72 = cr <8 3 ,73 = 
(7 2 (8 <7 3 ,75 = (7i <g> cr 3 . The s-wave superconducting OP 
and the number operator for the surface states of TIs as- 
sume identical form as for graphene, only without the <7o 
in the first block. Due to Nambu's particle-hole doubling, 
the true number of fermionic degrees of freedom on the 
surface of TI is a quarter of the one in graphene, how- 
ever. In this special case the critical theory acquires the 
supersymmetry, and the one-loop e-expansion is known 
to be exacts 

Next we wish to study the quantum phase transition 
from the semimetallic into the s-wave superconducting 
phase. Since we want to formulate an e(= 4 — d)- 
expansion near four space-time dimensions, we need 
first to define a spinor basis in which the theory can 
be formally extended from the physical three, to four 
space-time dimensions. We therefore rotate the spinor 
<J> — > U^f?, where, in graphene, U — exp [i f <7 ® 73] . Af- 
ter this unitary transformation, the s-wave OP reads 

$(cc) = (* 1 '(x)(cro<87ocos(p+cro ( 8*7o75sinv3)vl'(a;)), (2) 

while leaving the relativistically invariant free Dirac La- 
grangian, Lp unchanged. The number operator is then 

N = ctq ® 75. Similarly, in TI, the analogous trans- 
formation is performed by choosing the simpler U = 
exp [if 73]. 

For generality, we consider the U(f) gauge theory for 
Nf flavors of four-component Dirac fermions coupled to 



the bosonic OP with TVj, complex components via the 
Yukawa coupling in the presence of a fluctuating gauge 
field, with the complete Lagrangian L = Lf + Lb + Lbf + 
Lem- The coupling of the fermions to the U(l) gauge 
field reads 

L f = ^{x) lli {d ll -ie^A ll )if!{x). (3) 

The matrix 75 appearing in the minimal coupling is then 
the number operator, and e is the U(V) charge. The cou- 
pling of the OP to the masslcss fermions has the Yukawa 
form 

L bf = . 9 [(Re$)*# + (Imfc)*^*]. (4) 

On the other hand, the dynamics of the OP coupled to 
the U(l) gauge field can be described by the standard 
Ginzburg-Landau Lagrangian 

L b = + 2teA f ,)<f>\ 2 + m 2 |$| 2 + ^|$| 4 , (5) 

where m 2 is the tuning parameter of the transition. The 
U(l) gauge field is described by the usual Maxwell La- 
grangian 

Lem = -^FfivF^v, (6) 

with F^ v = d^A u — d u A^. We will use the transverse 
(Landau) gauge d^A^ = 0, in which the general gauge- 
invariant OP becomes locals The above theory is con- 
structed to be invariant under the following local U(l) 
gauge transformation <3> -> e le6 ' 75 ^', <I> — > e _2je6l $, A^ -> 

a^ + d^e. 

Without the gauge field (e = 0), the above field theory 
for Nf = 2 and Nf, = 1 also governs the critical behavior 
of graphene close to the transition to spin-singlet Kekule 
insulator. However, in the latter case, the Dirac spinor 
needs to be redefined as "J" = [^ + , ^_] T , where 'J a = 
[u CT (K + q), v a (K + q), u a (-K + q), v„{-K + q)], with 
a = ± as the spin projections, and with the frequency 
label suppressed.— We have set the Fermi velocity and 
the velocity of the bosonic excitations to be equal, since a 
weak anisotropy in the velocities is irrelevant^ The local 
(e 7^ 0) U(l) gauge theory describes the ultimate critical 
behavior at the superconducting transition, at which all 
the velocities in the theory are equal to the velocity of 
light. In graphene and on a surface of TIs, however, such 
a fixed point is experimentally inaccessible, since the bare 
Fermi velocities are ~ 10 6 m/s. 

Next, we proceed with the analysis of the £7(1)- 
symmetric Yukawa field theory. The couplings A, e and 
g are all dimensionless in (3+1) space-time dimensions, 
suggesting the e-expansion about d + 1 = 4 as a tool 
of choice for a study of the quantum critical behavior. 
Define then the action S rc „ = J dr J d d xL rcn , where the 
renormalized Lagrangian is 

7 rcn = ZvLf + Z*\(d f , + 2teA fl )<f>\ 2 + Z m m 2 \<f>\ 2 
+ Z x ^\<f>\ 4 + Z g L bf + Z A L E M- (7) 



The computation of the self-energy diagrams for the 
fermions, the order-parameter, and the gauge field using 
minimal-subtraction scheme then yields the renormaliza- 
tion constants to the one-loop order 



1 ,1 



Z* = l- -g 2 -, Z$ = l- g'Nf- + \2e z 



N h 



(8) 



(9) 



is, however, a relevant coupling at this critical point, 
(iii) the bicritical point in the e 2 = 0-plane: (g 2 ,A) = 
(e/X, (a — b)e), located in the unphysical region (A < 0) 
of the <I> 4 -interaction. Therefore, our one-loop results 
suggest that the semimetal-superconducting transition 
is of the second order in the absence of the fluctuating 
gauge field. The result is qualitatively similar to the in- 
sulating Ising and Hciscnberg universality classes*^ The 
correlation-length exponent (y) can readily determined 
from Eqs. © and iflj]). yielding 



where e = 4 — d and the dimensionless couplings Q 



{e 2 ,g 2 ,X} are rescaled as QSd/(2ir) 



), withS d = 



27r rf / 2 /r(d/2). (See Supplementary Material). The com- 
putation of the vertex diagrams to the same order gives 
the following renormalization conditions for the coupling 
constants: 



ry ryl/2 _ £ /2 



9 1 

3e 2 5 - 



9, 



(10) 



ZlXo^' - X 2 (N b +4)- - 96e 4 i + 2g i N f - = A. (11) 

The renormalization of the tuning parameter (to 2 ) can 
be extracted from the self-energy diagrams of the OP, 
leading to 



Z$TOg/x e 



X(N b 



l)-m 2 
e 



(12) 



Here, the couplings with subscript "0" are the bare cou- 
plings, and the ones without the subscript are the renor- 
malized couplings, and (i is the renormalization scale. 
Dimensional regularization explicitly preserves gauge in- 
variance of the theory implying /i^Z^e 2 , = e 2 , to any 
orderi 18 i 22 In conjunction with this identity and Eq. ©, 
one can write the (ultraviolet) beta-function of the charge 



*L—* + ^s, + »y. (i3, 



The renormalization group flow of the remaining two cou- 
plings can be obtained from Eqs. (|TU|) and pT|) 



-eg 2 + (N f + l)g A - I8e 2 g 2 , (14) 
-eA + 2N f g 2 (X - g 2 ) - 24e 2 (A - 4e 2 ) 



(N b + 4)A 2 



(15) 



The above beta-functions, besides the trivial, yield 
the following neutral (e = 0) fixed points: (i) the 
Wilson-Fisher fixed point at g 2 = and A» = 
e/(N b + 4); (ii) the neutral Gross-Neveu fixed point: 
( g 2 A,) = (e/X, (a + b )e) , where a = (1 - N f )/2XW, b = 
y/(N f - l) 2 + 8N f W/2XW, with X = N f + 1,W = 
Ni, + 4. Ignoring the gauge coupling e for the mo- 
ment, this fixed point is critical, and it controls the 
transition towards the s-wave superconducting state, or 
into the spin-singlet Kekule insulator. Weak charge e 2 



v=\ + \{N b + l)X^-iel + ^gl 



(16) 



with e 2 = and (g 2 , A* ) corresponding to the neutral 
Gross-Neveu critical point. Since the Lorentz-symmetry 
breaking perturbations are irrelevant near the critical 
points the dynamical critical exponent is z = 1, and 
the Fermi velocity (vf) is non-critical. Near the neutral 
critical point both the order-parameter and the fermion 
fields acquire a non-trivial anomalous dimension which, 
respectively, reads 



l2el)e + <D(e 2 ), r, f 



^e + 0(e 2 ). (17) 



The residue of the quasiparticle pole of the fermions 
Z ~ m ZU71f ~ TO £ / 4;i<: , and they cease to exist as sharp 
excitations at the critical point. Moreover, as the system 
approaches the critical point from the superconducting 
side both the mass of the superconducting OP and the 
fermion mass vanish with the universal ratio 



2 A* 
5* 



(18) 



In order to extract the critical exponents and amplitudes 
for graphene one needs to substitute Nf = 2, N b = 1, 
while for surface states of TIs one should use Nf = 1/2, 
N b = 1. In the latter case, we obtain Vh = f]f = e/3 and 
^ = l/2 + e/4in agreement with Ref. [20|. 

In the fully gauged theory with e ^ 0, the 
charged Wilson-Fisher fixed points are at (e 2 ,g 2 , Xf) = 



„ i8+y±V(is+y) 

' U ' 2YW 



2 -216W 



e, where Y = Nf + N b . 



On the other hand, previously discussed neutral Gross- 
Neveu fixed point is unstable in the charge direction, and 
a pair of charged fixed points is located at 



(e*, <?*, Xf) 



3 27 
4F" " 



2Y Ax ± ^A 2 + A 2 



2XY 



2X 2 Y 2 W 



(19) 

where A i = XY [XY + 18 X - Nf (27 + 2Y)] , A 2 = 
-4WX 2 Y 2 [54X 2 -27V/(13.5 + y) 2 ]. However, only 
the fixed point with A+ > is stable in the critical plane 
(m 2 = 0). This fixed-point is therefore critical, and con- 
trols the behavior in the vicinity of the quantum phase 
transition in the full Lorentz-invariant U(l) gauge the- 
ory. Furthermore, for any physical number of flavors this 



critical point describes the second-order phase transition, 
since the quantity A 2 + A2 is positive for any Nf > 0.142. 
However, since all the velocities in this theory are set to 
be equal to the velocity of light, this critical point may 
be reached only in the deep infrared regime^ The other 
fixed-point at (e 2 ,g*,A~) lies in the unphysical region 
(A < 0) of the $ 4 -interaction for any Nf. On the other 
hand, when Nf = 0, we obtain the standard one-loop 
result for the critical number of the complex components 
of the OP above which the superconductor transition 
is of the second order, Ng rq,t ~ 182.952, and the tran- 
sition is controlled by the charged Wilson-Fisher fixed 
point, with g 2 = 0i 18 i 22 It is worth observing that with- 
out the Yukawa interaction, N£ rit reduces to 3.47 if one 
takes only into account the coupling of fluctuating gauge 
fields with the massless Dirac fermions<2i The Yukawa 
coupling therefore appears to be crucial for the stabiliza- 
tion of criticality in the theory, and the suppression of the 
possible discontinuous transition, which occurs in related 
theories J£ 

The superconducting coherence length (£) diverges as 
£ ~ m~ v , and the correlation length exponent [y) can 
be computed readily from Eq. (|16[) . The boson and 
the fermion fields in the vicinity of this charged critical 
point acquire anomalous dimensions, which can be found 
from Eq. (TTT)) . One can also compute the flow for the 
Ginzburg-Landau parameter k 2 = A/(2e 2 ) characterizing 
the transition, reads 



3 Y 



12 \k 



2N f 



(20) 



At the charged Gross- Neveu critical point, this flow equa- 
tion has fixed points at k 2 _ < and k 2 , > 0, for ar- 
bitrary Ni, and Nf > 0.142. The residue of the quasi- 
particle pole vanishes at the charged critical point as 

The optical conductivity in the entire semi-metallic 
phase remains constant and universal, while it becomes 
infinite (zero) in the superconducting (Kekule) phase. 
Right at the quantum critical point it is also expected 
to be universal, but different from the one in the semi- 
metallic phased The universal conductivity at the Gross- 
Neveu (neutral or charged) critical point is also expected 
to be different from one found in a pure bosonic theory— 
The computation of its value is an interesting problem 
left for future research. 

To summarize, by employing a [/(l)-symmctric 
Gross-Neveu- Yukawa theory, we here studied the zero- 
temperature scmimetal-superconductor (Kekule insula- 
tor) transition in graphene and surface of TIs, and 
showed that it is continuous for any number of Dirac fla- 
vors. The full U(l) gauge theory exhibits a charged crit- 
ical point also for an arbitrary number of Dirac flavors, 
and may be relevant for the semimetal-supcrconducting 
transition in the deep infrared regime. 
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We here present details of the renomalization group calculation to the one-loop order using the Euclidean partition 
function^ For the ease of calculation, we here rewrite the coupling of the order parameter to the massless Dirac 
fcrmions, shown in Eq. (4) of the main text as 

L bf = .g[(Rc$)§* + (Im$)*i75*] = + (21) 

where the projectors are P± = |(1 ± 75). 



I. COMPUTATION OF BOSONIC SELF-ENERGY 



First we consider the renormalization of the self energy of the order parameter field, arising from the diagrams 
shown in Fig. [TJ The one loop correction of the bosonic self energy due to its coupling with the gauge field (see Fig. 

< 2ie > 2 /!?? ( v - v -f) p* + pm* + >). (k + +nf -*? a - h\ , (22) 



QJa)) reads 

(la) = - 



where 



and 



j _ f d d P {2k + P f _ r 1 r d d P (2-x) 2 fc 2 + P 2 

1 J (2tt)«» p 2 p+p) 2 +m 2 ] J * J {2vY \p* +x (l- x )k*+xm?f [ ' 



, 2 

U2 



T _ f d d p p,{2k + p) ll p v {2k+p) u f 1 f d d p ((p-xk) 2 + 2k-p-2xki / 

H J (2tt)<* p±[(k+p)*+m*\ J aXZ[l X) J (2tt)<* tf+ x (l-x)k*+m*] a [ ' 

after setting p + X k — > p in both the integrals. Performing some standard integrals we arrive at 

I 1 = ^(2k 2 -m 2 )+0{l), / 2 = ^(_fc 2 - m 2 )+0(l). (25) 

Therefore the correction to the bosonic self energy due to its coupling with the gauge field is 

(la) = 12e 2 fc 2 ^+0(l). (26) 

Besides the above renormalization, the bosonic self energy acquires additional correction from the interaction vertex 
(A|$| 4 ), as shown in Fig. [IJb). This contribution is 

A P d d a 1 AT 

(16) = - (2N b + 2) / — L - T — - = A (Nb + 1) — - m 2 - + 0(1). (27) 
2 J {2TT) a q z + m z e 

Apart from the above two contributions, additional renormalization of the bosonic self energy comes from its coupling 
with the massless Dirac fermions, as shown in Fig. [He), yielding 



where d = 7 M g M . These three diagrams together give the following renormaization conditions 



+g 2 N f ^ - 12e 2 ^ = 1 (29) 
Z$m^- £ - A(7Vf, + l)m 2 i =m 2 , (30) 



leading to Z$ in Eq. (8) and mass renormalization in Eq. (12) in the main text. 
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FIG. 1: Diagrams contributing to the self energy corrections of the bosonic order parameter field. 
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FIG. 2: Diagrams contributing to the self energy corrections of the gauge field. 



II. SELF ENERGY CORRECTION OF THE GAUGE FIELDS 

Next wc present the computation of the self energy correction of the gauge fields. It arises from the diagrams shown 
in Fig. [3] Contribution from the diagram (a) is 



(2a) = N b {2ef 



2 f d d q (k + 2q)^(k + 2q) L 



= 4 N h e 2 



dx 



d d q | q 2 5 p „ + (1 - 2x) 2 k p k v 



(27r) d [q 2 + m 2 ] [(q + k) 2 + m 2 ] ' J ™ J (2tt)«* [ q * + x (\ - x )k* + m 2 ] 2 

once we take q + xk — > q. Performing some standard integrals gives 



(2a) 



2 Nd ,2/r k p k v 
N b e — k [ dp,, 



8A b e 2 ^ o^m 2 + 0(1). 



(31) 



(32) 



Upon taking the contribution from the diagram in Fig. [21c), the term proportional to m 2 in the last expression exactly 
cancels out and the gauge field remains transverse. The non-trivial coupling of the gauge field with Dirac fermions 
also gives a correction to the gauge field propagator, as shown in Fig.[^b). Contribution from this diagram reads as 



(26) = e 2 Tr 



d d q (f + <f) 



where 



(27t) d q 2 (k + q) 2 

d d q q p (k + q) 



e 2 Tr [7^ j u lo] I pa = e 2 AN f (6^5^ - 5^5^ + 5^5^) I p<T , (33) 
d d q q P q a - x(l - x)k p k a 



= / dx 



(2n) d q 2 (k + q) 2 J J (2i:) d [ q 2 + x{1 _ x)k 2f ' 
once we take q + xk — > q. A tedious, but otherwise straightforward calculation yields 



0(1) 



Therefore the correction to the gauge propagator from the diagram in Fig. EJb) is 



(26) = ^ A /e 2 k 2 



N d k- 



0(1). 



Collecting the contributions from these three diagrams gives the renormalization condition for Z a as 

Za + e 2 \ (N b + N f ) ^ = 1, 
leading to the announced result in the main text, Eq. (9). 
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FIG. 3: Diagrams contributing to the self energy corrections of the Dirac fermions. 
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FIG. 4: Renormalization of the fermion-boson vertex. 



III. CORRECTION TO THE SELF ENERGY OF DIRAC FERMION 



Let us now compute the fermionic self-energy corrections to the one loop order. Fermions coupling with the gauge 
field, shown in Fig. O^a), leads to the following correction to its self-energy 



(3a) 



d d P i{f + $) 1 / p„p v 



(2 - d) 7p 
d d p 



d d p (p + k) p I" d d p p fJ ,(k+p) p p 1/ 

{2ir) d p 2 {k+p) 2 ^ lpl "J (2n) d p*(k+pf 



{2tt)° 



[p 2 + x(l — x)k 2 



[p 2 + x(l — x)k 2 



(38) 
(39) 



where p + xk — > p. After performing the integrals we find a trivial contribution 



(3a) 



-i e 



(2 - d) - - (2 - 2d) 



(4tt) : 



+ 0(1) = + 0(1), 



(41) 



as d — > 4. However, the fermionic self energy acquires nontrivial correction due to its coupling with the order-parameter 
field, shown in Fig.^b), yielding 



(36) = g 2 



d d k 
(2+F 



i 



(fc — p) 2 + m 2 



1 



9 \Y e ) ifN d m- e + 0{l). 



(42) 



Therefore 



= 1 - — Nd m" 



(43) 



as shown in Eq. (8) in the main part of the paper. 



IV. RENORMALIZATION OF BOSON-FERMION VERTEX 



To the one loop order the boson-fermion vertex also gets renormalizcd. Its renormalization comes from the diagrams 
shown in Fig [4] Let us first consider the contribution of the diagram (a), renormalizing P + vertex 



(4a) 



(«0 2 
-P+e 2 



d d q 1 



(27r) d q 2 
d d q 
(2n) d 



-n [Sau- 



p 



(k - q) 2 ^ (k + k> - qf 



l5lv 



where 



d d q (k-q) p (k + k' -q) a 



and 



(27r) d q 2 {k - q) 2 (k + k' - qf 



d d q {q-k) p (q-p) a q ll q v 



= / Zydy 



Iplplalplpa ~ lplplalvlp.vpa 

d d q [q + (l- v)p - (1 - xy)k) p (q + xyk - yp) 



(44) 
(45) 



dx 



(2tt) c 



((1 — y)p + xyk) + xyk 2 + (1 — y)p 2 



(46) 



(27T) d q A {q-k) 2 {q-pf 



= 6 



y dy I (1 — x)dx 



d d q 
(2n) d 



Qp(q,P,k,x,y)Qv(q,p,k,x,y) [q+ (1 - y)p + (xy - l)k] [q - yp + xyk] c 



((1 — y)p + xyk) + (1 — y)p 2 + xyk 2 



(47) 



after shifting the momentum as q — (1 — y)p — xyk q in I p(T , Ipvpa, and Q(q,p, k, x, y) = q + (1 — j/)p + zy/c. For our 
renormalization group calculation, since we use the minimal-subtraction scheme, we only need to keep the divergent 
pieces from these two integrals, which read 



jdiv 
pa 



fdiv 



N d S pa — fc" e , 

6 / y 2 (l — x)dxdy 



(48) 



d d q 
{2it) d 



qpQvqpQa 



q 2 + (1 — y)p 2 + xyk 2 — ((1 — y)p + xyk)' 



N d 
24 e 



Hence the renormalization of the boson-fcrmion vertex is 



(4a) 



-P+e 2 



IplplalpSpa^ ~ Iplplalv {SpvSpa + ^pp^va + Spa$pv) ^ 



■0(1) 



-P. 



3e 2 



(49) 



JV d + 0(l). (50) 



The renormalization of the other boson-fcrmion vertex, proportional to P_ , comes from the diagram (b) of Fig. |4l 
which gives an identical contribution as in Eq. |50j . In principle, the boson-fermion vertex can also be renormalized 
from the diagram Fig. H](c) and (d). They respectively renormalizes the P+ and P_ vertex. However, it is sufficient 
to consider one of them, for example (c). Its contribution reads as 



(4c) = -(ie)(2ie) 



= -2P+e' 



d d p 1 
(2n) d p 2 



Pp Pv 

J) 2 



(2q + p), 



(p + q) 2 + m 2 (q+p + k)'- 



P+ 



757^ 



d d p 
(2^p 



(51) 



where 



I is fit — 



d d p 



{2q + p) v (q + k') a 
(27r) d p 2 [(p+q) 2 +m 2 } (p + k') 2 



2y dy 



dx 



d d p (p + (1 + y)q - xyk) v (p - (1 - y)q + (1 - xy)k') a 



(2tt) c 



b 2 + A]- : 



(52) 
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FIG. 5: Renormalization of the interaction vertex of boson. 



where, k' = k + q and A = xyk' 2 + (1 — y) [q 2 + m 2 ) — [(1 — y)q — xyk'] 2 . We shift the momentum p+ (l~y)q+xyk' 
p, while writing the last equation. The other integral I^^a reads as 

d d p p^(2q+p) (p + k') a 



f d a p 
W« - J (27r)d p 



4 l(p + q) 2 +m 2 } (q + k') 2 

- 1 2 , r 1 ,, H fl (x,y,q,p,k , )H^x,y,q,p,k')(p+{l + y)q-xyk') {p-(l-y)q-{xy-l)k') c 
o y ay (1 — x) ax 



[p 2 + A] 4 



(53) 



after shifting the momentum p as before, where H(x, y,q,p,k') = p — (1 — y)q — xyk' . We are however interested only 
in the divergent pieces of these two integrals which read as 



I, 



rdiv 



^k-s 

4 e 



N d 
24 e 



(54) 
(55) 



Hence the renormalization of the P+ vertex from the diagram (c) in Fig. 0] reads as 



(4c) = -2 P+ e 2 



lortv S aL , 7 a 7„ 



24 



— + 0(1) = + 0(1), 



(56) 



as <i — > 4. Hence this diagram provides trivial renormalization of the boson- fermion vertex (P+). Similarly, one can 
show that renormalization of the other boson- fermion vertex (P_ ) arising from the diagram Fig. 0] (d) is also trivial. 
Renormalization of the Yukawa vertex can also arise from boson- fermion interaction, as shown in Fig. [4] (e), (f), (g), 
(h). The combined contribution of these diagram is however 



(4e) + (4/) + (4 ff ) + (4/i) = g 



2 / d d q 



P+ 



P+ 



(27r) d 



(k+p-qf (k-qf 



P_ 



P_ 



P+ 



P4 



(k + p + qf i+ (fc + g) 2 + 



(k-p-q) {k-qf 



P_ 



P_ 



*(jS-jf + f) + 



P4 



(fc-p + g) 2 ' (fc + g) 2 + 



g 2 + to 2 



= 0. 



(57) 



These contributions yield the announced renormalization condition of the Yukawa coupling (g) in the main text, 
namely 



Z4, gofJ, + 3e g - = g. 



(58) 



V. RENORMALIZATION OF |$| 4 VETREX 

|4 



Finally we consider the renormalization of the |$| vertex. Renormalization of this vertex from the Yukawa coupling, 
shown in Fig. E](a), reads 



(27r) d 



q 2 (pi + q) 2 (pi+q- Pi) 2 {jpi + q - Ps - Pif 



The rcnormalization of this vertex arising from the diagrams in Fig. [5fb) is 



(5b) ^X 2 (N b + 4) ^ + 0(1) 



(60) 



Moreover, the coupling of the order parameter field with the gauge field also provides a renormalization of the |<I>| 4 
vertex, as shown in Fig. [SJc). Its contribution reads 



(5c) 



32 e 4 



d d q 1 



(27r) d g 4 



Suu - 



<v - 



= -32 e 



4 37V d 



(61) 



Upon collecting contributions from all the diagrams in Fig. [SJ one gets the renormalization condition of the coupling 
constant A, as in Eq. (11) in the main part of the paper. 
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